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Abstract. The paper is concerned with the equilibrium distributions of 
continuous-time density dependent Markov processes on the integers. These 
distributions are known typically to be approximately normal, with 0{1/ ^/n) 
error as measured in Kolmogorov distance. Here, an approximation in the 
much stronger total variation norm is established, without any loss in the 
asymptotic order of accuracy; the approximating distribution is a translated 
Poisson distribution having the same variance and (almost) the same mean. 
Our arguments are based on the Stein-Chen method and Dynkin's formula. 



1. Introduction 

Density dependent Markov population processes, in which the transition rates 
depend on the density of individuals in the population, have proved widely 
useful as models in the social and life sciences: see, for example, the monograph 
of Kurtz (1981), in which approximations in terms of diffusions are extensively 
discussed, in the limit as the typical population size n tends to infinity. Here, we 
are interested in the behavior at equilibrium. Our starting point is the paper of 
Barbour (1980), in which conditions are given for the existence of an equilibrium 
distribution concentrated close to the deterministic equilibrium, together with 
a bound of order 0(1/ ^/n) on the Kolmogorov distance between the equilibrium 
distribution and a suitable normal distribution. We now show that this normal 
approximation can be substantially strengthened. Using a delicate argument 
based on the Stein-Chen method, we are able to establish an approximation 
in total variation in terms of a translated Poisson distribution. What is more, 
our error bounds with respect to this much stronger metric, and under weaker 
assumptions than those previously considered, are still of ideal order 0{1/ \/n). 

The first step in the argument is to establish the existence of an equilibrium 
distribution under suitable conditions, and to show that it is appropriately con- 
centrated around the 'deterministic' equilibrium, defined to be the stationary 
point of an associated system of differential equations which describe the aver- 
age drift of the process in the limit as n — s> oo; this is accomplished in Section [21 
The closeness of this distribution to our approximation is then established in 
Section HI by showing that Dynkin's formula, applied in equilibrium, yields an 
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equation not far removed from the Stein equation for a centred Poisson distri- 
bution, enabling ideas related to Stein's method to be brought into play. An 
important element in obtaining an approximation in total variation is to show a 
priori that the equilibrium distribution is sufficiently smooth, in the sense that 
translating it by a single unit changes the distribution only by order 0(1 /y/n) 
in total variation: see, for example, RoUin (2005). The corresponding argument 
is to be found in Section [3l We illustrate the results by applying them to a 
birth, death and immigration process, with births occurring in groups. 

1.1. Basic approach. We start by defining our density dependent sequence of 
Markov processes. For each n G N, let t > 0, be an irreducible continuous 

time pure jump Markov process taking values in Z, with transition rates given 
by 

i ^ i + j at rate nXj ^— j , -i G Z, j G Z \ {0}, 
where the Aj(-) are prescribed functions on M; we set 

Zn{t) := n-'Zn{t), t>0. 
We then define an 'average growth rate' of the process Zn a.t z & n~^Z by 

F{z) := ^ jX,{z), 

jez\{o} 

and a 'quadratic variation' function by n~^a'^{z), where 

j&\{0} 

assumed to be finite for all 2 G M. 

The 'law of large numbers' approximation shows that, for large n, the time 
dependent development of the process Zn runs close to the solution of the dif- 
ferential equation system z = F{z), with the same initial condition, and that 
there is a approximately diffusive behaviour on a scale n~^^^ about this path 
(Kurtz 1970, 71). If F has a single zero at a point c, and is such that c is glob- 
ally attracting for the differential equation system, then Z„ has an equilibrium 
distribution n„ that is approximately normal, and puts mass on a scale n^/^ 
around nc (Barbour 1980). The corresponding asymptotic variance is given by 
n^/'^Vc with Vc := ^2F'\c) ' Pi^o^ided that F'{c) < 0, and the error of the approx- 
imation in Kolmogorov distance is of ideal order 0(?t.~^/^) if only finitely many 
of the functions Aj are non-zero. 

In this paper, we strengthen this result, by proving an accurate approxima- 
tion to the equilibrium distribution using another distribution on the integers. 
Under assumptions similar to those needed for the previous normal approx- 
imation, we prove that the distance in total variation between the centred 
equilibrium distribution n„ — [ncj and the centred Poisson distribution 

Y'o{nVc) ■■= Po(m;c) * 

is of order 0{n~^/'^): here and subsequently, 5r denotes the point mass on r, 
and * denotes convolution. If infinitely many of the \j are allowed to be non- 
zero, but satisfy the analogue of a (2 -|- a)'th moment condition, for some 
< a < 1, we prove that the error is of order 0(r2~"/^). 
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The proof of our approximation runs as follows. The infinitesimal generator 
An of Zn, acting on a function h, is given by 

iAnh){i) := nXj(^-^[h{i + j)-h{i)], i e Z. 

jez\{o} 

In equilibrium, under appropriate assumptions on h, Dynkin's formula implies 
that 

(1.1) E{Anh){Zn) = 0. 

The following lemma, whose proof we omit, expresses Anh in an alternative 
form. 

Lemma 1.1. Suppose that X]jez\{o} ^ ^ /'^^ -2 G M. Then, for any 

function /i: Z — M with bounded differences, we have 

(1.2) {Anh){i) = VQ) v^7h(0+^^(^)^7h(0 + K(^7,0, 
where XJ f{i) := f{i) — f{i — 1) and Qhif) '■= \/h{i + 1) and, for any i G Z, 

(1.3) := -^fQ) v^/.(0 + Xl^j(^'')^^j(~) ~Z1^J(^'^)^^-J(~ 



i>2 j>2 



with 



;i.4) 2a,{g,i) := - 1) g{i) + 2^,^ \7 9{i + J - k) 

k=l 

3 



(1-5) = 2^Q^ v'^7h(^ + j-/^ + i); 

26,((7,z) := j{j-l)yg{i)-2Y^k\jg{i-j + k) 

k=l 

= 2X^(2) V'^7h(^-J + ^)- 

fe=2 ^ ^ 

Writing (11. II) using the result of Lemma 11.11 leads to the required approxi- 
mation, as follows. In equilibrium, Zn/n is close to c, as is shown in the next 
section, and so the main part of (11. 2p is close to 

-F\c) 1^^;^ V 9H{i) -{^- nc)g,{i) 

because F{c) = 0. Here, the term in braces is very close to the Stein operator 
for the centred Poisson distribution P{nVc) with = applied to the 

function g^. see Rollin (2005). Indeed, for any f > and B C Z^,, where 
Z„ := {Z G Z, / > — \v\ }, one can write 

(1.6) Ml) - Vo{v){B} = V ^g{l + l)- Igil) + {v)g{l), I e Z„, 
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for a function g = Qv b satisfying 
(1.7) 

sup |^7(/+1)| < minjl,^); sup |v^7(/+l)l < -; =0, /< 

l>-[v\ ^ \/vi l>-\v\ V 

where (x) := x — \_x\ denotes the fractional part of x; note also, from (11. 6p 
and ([LTD, that 

(1.8) sup|/(7(/)| <3. 

Replacing / in (11.61) by an integer valued random variable W then shows that, 
for any B C Z^, 

\¥[W G B]-F^{v){B}\ 

(1.9) < snp \K{v \/ g{W +1) -Wg{W) + {v)g{W)}\ + ¥[W < -[v\], 

g&Gv 

where denotes the set of functions Z — M satisfying (ll.7p and (11.81) . 
Hence, replacing W by Z„ and v by nvc in (11.91) . and comparing the expectation 
with (II. ip expressed using Lemma II. H the required approximation in total 
variation can be deduced; for this part of the argument, we need in particular 
to show that, in equilibrium, 

(1.10) |E{v^?(^„ + l)- V^7(^n)}| = |E{v'^^(^n+l)}| = 0(n-3/2), 

and also that W\En{g, Zn)\ = 0{n~"^'^) for any g E Qnvc- The bound (ll.lOp 
follows from Corollary 13.31 in Section [3], and the latter estimate, which also 
uses (1 1.1 op . is the substance of Section HI 

1.2. Assumptions. We make the following assumptions on the functions Xj. 
The first ensures that the deterministic differential equations have a unique 
equilibrium, which is sufficiently strongly attracting. 

Al: There exists a unique c satisfying F(c) = 0; furthermore, F'{c) < and, 
for any t] > 0, fir^ := inf|^_c|>^ 1-^(2^)1 > 0. 

The next assumption controls the global behaviour of the transition func- 
tions Xj. 

A2: (a) For each j G Z \ {0, }, there exists Cj > such that 

(1.11) Xj{z) <Cj{l + \z-c\), zeR, 
where the Cj are such that, for some < a < 1, 

j€Z\{0} 

(b) For some A° > , 

Xi{z) > 2A°, zeR. 

The moment condition on the cj in Assumption A2 (a) plays the same role as 
the analogous moment condition in the Lyapounov central limit theorem. Un- 
der this assumption, the ideal rate of convergence in the usual central limit 
approximation is the rate 0(n~"/^) that we establish for our total variation ap- 
proximation. Assumption A2 (b) is important for establishing the smoothness 
of the equilibrium distribution n„. If, for instance, all jump sizes were multiples 
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of 2, the approximation that we are concerned with would not be accurate in 
total variation. 

We also require some assumptions concerning the local properties of the func- 
tions Xj near c. 

A3: (a) There exist s > and < 5 < 1 and a set J C Z \ {0} such that 

inf Xj{z) > eXj{c) > 0, j E J; 

\z—c\<S 

Xj{z) = for all \z - c\ < 6, j ^ J- 
(b) For each j G J, Xj is of class C"^ on \z — c\ < 6. 
Assumptions A2 (a) and A3 imply in particular that the series Yljez\{o} 
and '^j^z\{o}j^-^ji^) uniformly convergent on \z — c\ < 6, and that their 
sums, F and respectively, are continuous there. They also imply that 

^ \j\nXj{i/n) = 0{\i\), \i\ oo, 

i6Z\{0} 

so that the process Z„ is a.s. non-explosive, in view of Hamza and Klebaner (1995, 
Corollary 2.1). 

The remaining assumptions control the derivatives of the functions Xj near c. 
A4: For S as in A2, 

T WKWs 

Li := sup < oo, 

where := sup|^_^|<5 \f{z)\. 

This assumption implies in particular, in view of Assumptions A2-A3, that 
the series X]jGZ\{o} ^'^i(^) Sjez\{o} ^^'^j(^) uniformly convergent on 
|z — c| < 6, that their sums are F' and (cx^)' respectively, and that F and o"^ 
are of class on \z — c\ < 6. 



A5: For 6 as in A2, 



I 

\Xj\\5 



This assumption implies, in view of A2-A3, that the series J2jez\{o} 
uniformly convergent on \z — c\ < 6, its sum is F", and F is of class on 
\z — c\ < 6. 

Our arguments make frequent use of the following theorem, which is a re- 
statement in our setting of Hamza and Klebaner (1995, Theorem 3.2), and 
justifies f ll.ll) . 

Theorem 1.2. Suppose that Zn is non-explosive. Let h be a function satisfying 
(1.12) 

{\An\h){t) := < Cn4iy\h{t)\), H^oo, 

jez\{o} ^ 

for some Cn,h < oo. Then, if h{Zn{0)) is integrable, so is h{Zn(t)) for any t > 0; 
moreover, 

/l(Z„(t))-/l(Z„,(0))- / {Anh){Zn{s)))ds 
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is a martingale, and Dynkin's formula holds: 



(1.13) E[h{Zn{t)) - h{Zn{0))] = [ E{Anh){Zn{s))ds. 

Jo 

2. Existence of the equilibrium distribution 

In this section, we prove that Z„ has an equihbrium distribution which is 
suitably concentrated in the neighbourhood of nc. 

Theorem 2.1. Under Assumptions AI-A4, for all n large enough, Zn has an 
equilibrium distribution Un, and 

EnJI^n, - c| ■ a(|^„ - c| > 6)} = 0{n-') 

En„{(^„ - cf ■ a(|2„ - c| < 6)} = O(n-i), 

for 5 as in Assumption A3: here, as before, Zn '■= n~^Zn. 

Proof. The argument is based on suitable choices of Lyapounov functions. Con- 
sider the twice continuously differentiable function M — > IR^ defined by 
V{z) := \z — cp"*"", for the a in Assumption A2(a). Since V{c) = and 
V{z) > for any z ^ c, and because 

(2.2) F{z)V'{z) = -|F(2)|(2 + a)|^-c|^+" < for any ^ ^ c, 

while F{c)V'{c) = 0, we conclude that \^ is a Lyapounov function guaranteeing 
the asymptotic stability of the constant solution c of the equation x = F[x). 
We now use it to show the existence of n„. 

Lemma 2.2. Under the assumptions of Theorem \ 2.1\ the function hy{i) := 

^(n) ~ In ^ c|^^" fulfils the conditions of Theorem with respect to the 
initial distribution 61, the point mass at I, for any I G Z. 

Proof. Checking (11.121) . we use Taylor approximation and Assumption A2 (a) 
to give 

i\An\hvm < (2 + a)|z-c|^+" Yl IjI c,(l + k-c|) 

jez\{0} 

(2.3) ^ (2-H„)(l.„)|.-c|- ^ /c,(l + |.-c|) 

jez\{o} 

(2.4) T, brc,(l + |.-c|), 

iez\{0} 

where we write z := i/n. For l-z — c| < 6 < 1, the estimate in (12.31) is uniformly 
bounded by 

C„. := 2(2 + a){$:bte+(i±^$:fc, + ii±^5:i,rc,} < 

j j j 

because of Assumption A2 (a); for |;z — c| > 5, we have the bound 

as required. □ 



00, 
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The above lemma allows us to apply Dynkin's formula to the function hy- 
Using Taylor approximation as for (12. 3p . but now noting that the first order 
term 

^ Hz)jV\z) = F{z)V\z) 
ieA{o} 

can be evaluated using f l2.2p . it follows that 

(2.5) (AMW < -\F{z)\{2 + a)\z-c\^+'' + n"^C2 < n'^Ca 
on — c| < 5, for 

C, = (2 + a)(l + a){5^fc, + ^|jf+%.} < oo, 

j j 

where, once again, z := i/n. On |-2 — c| > 6 and under Assumption A2 (a), we 
have 



(Anhvm < -|F(2)|(2 + «)|2-c|i+° 

(1 + a) 



2n\F(z)\-\z-c\ 

I ^ ^1 I I jez\{o} 



— T E fc,{l + \z-c\) 



' ^ ^' ' ' iGZ\{0} 

(2.6) < -^^^^^±^\z - c\'+- < -f^s\z-c\'+-, 
as long as n is large enough that n6 > 1 and 

^ iez\{o} 
Dynkin's formula (11.131) then implies, for such n, that 

< EihviZnit)) = V{z)+ [ E,{Anhv){Zn{s))ds 

Jo 

< V{z)+ [ —¥i{\n~^Z4s)-c\<5)ds 
Jo n 

-lis I E,{|n-iZ„(s)-c|i+"-]l(|r2-iZ„(s)-c| >5)}cis, 
Jo 

for any t > and z G Z, where Pj and Ej denote probability and expectation 
conditional on Z„(0) = i. It now follows, for any y > S, that 

rP,(|n-iZ„(5)-c|>y)d. 

< ^/ E,{\n''Z^{s)-c\'+--l{\n-'Z^{s)-c\>y)}ds 
't Jo 

(2.7) < -V{z) + ^! ¥i(\n-^Zn{s)-c\<5)ds, 

t nt ./n 



Co 
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and, by letting t — oo, it follows that 
1 ft 

limsup- / Pj(|n ^ Zn{s) — c\ > y) ds < 

t— >oo t Jo 

This implies that a limiting equilibrium distribution n.„ for Zn exists, see for 
instance Ethier and Kurtz (1986, Theorem 9.3, Chapter 4), and that, writing 
Zn '■= n~^Zn, we have 

Pn„(kn-c| >!/) < ^\ 
for any y > S. Furthermore, 

PCX) 

EnA\zn - c\ ■ t{\zn - c\ > 6)} = FuAl^n - c\ > y) dy 

poo 

proving the first inequality in (12.11) . 

For the second inequality in (12. ip . we define a function V : M — M, which is of 
class C^(M), is bounded and has uniformly bounded first and second derivatives 
on M, fulfils the conditions of Theorem 11.21 and satisfies F{z)V'{z) = —{z — c)^ 
on l-z — c| < 6. 

In view of the latter property, we begin by letting v. [c — 6,c + 6] M+ be 
the function defined by 

v{z) := -^^dx, 

with v{c) = 0. Note that v is well defined, since F'{x) < on a small enough 
neighborhood of c, by Assumptions Al and A4, and that v{z) > for any z ^ c. 
Furthermore, in view of Assumptions Al and A4, 

^ - F{z) ^""^ ~ F^z) 

exist and are continuous on l^; — c| < 6, since |-F(2;)| > for z 7^ c, F{z) ~ 
F\c){z — c) for z c, and F' is continuous. In particular, we have 

(2.8) v'(c) = limf'(2;) = and v"(c) = \imv"(z) = — > 0. 

z-»c z-*c F'[c) 

Now define the function V to be identical with on 1^; — c| < 5, and continued 
m z < c — 5 and in 2; > c + 5 in such a way that the function is still C2 , and 
takes the same fixed value everywhere on l^; — c| > 25. Let 

C3 := maxjsup V^(2;), sup |V^'(2;)|, sup |y"(2;)|}. 

zeM 2eiR ^eM 

Lemma 2.3. Under the assumptions of Theorem \ 2.1[ the function hv{i) '■= 
^in) f'^^fi^'^ the conditions of Theorem M.^ with respect to the initial distribution 

n„." 
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Proof. Since hv{i) is bounded, it follows that Ku„\hv{Zn)\ < oo. |.4.n| /iy is 
also bounded, since, for \n~^i — c\ < 46, by Assumption A2 (a), 

{\An\hvm < C3 c,(l + 45), 

j£Z\{0} 

while, for \n~^i — c\ > 46, 

{\An\hvm < Cs Yl c,{l + \n-H-c\) 

j : \j+i—nc\<2n5 
i6Z\{0} ' ' i6Z\{0} 

□ 



We now apply Dynkin's formula to hy, obtaining 

■2 

= En„{(AM(^n)} < En„{F(2„)V^'(^„)+ Yl ^J^^^^i^ ^^j' 

jez\{o} ^ 

Hence it follows that 

EnA-F{zn)V'{zn)-l{\zn-c\ <6)} 

■ 2 

< En„{F{zn)V'{z^)-l{\z^-c\>6)+ J] A,(z„)l- C3}, 

iez\{o} ^ 

whence we obtain 

EuAiZn - CY ■ l{\Zn - C\ < 6)} 



2 

< EnA\F{zn)V'{z^)\-l{\z^-c\>6)} + CsEn„{ Y ^^^'^^^l) 

< Cs Y ('^\j\ + -)^J^TiAK-c\-li\z^-c\>6)} + -^ sup a'iz). 



j&Z\{0} 

Using the first inequahty in (12.11) and Assumptions A2 and A3, we conclude 
that 

En„{izn-cf-mzn-c\<6)} = Oin-'), 
proving the second inequality in (12. II) . □ 

Corollary 2.4. Under Assumptions AI-A4, 

EuAK-c\} = Oin-'/'). 

Proof. Using Holder's inequality, we obtain 

E{\Zr,-c\} 

= En„{|2;„ - c| ■ l{\zn -c\> 6)} + EuAl^n - c\ ■ l{\zn -c\< 6)} 
< E{\z^ - c\ ■ l{\zn -c\> 6)} + VEuAi^n " c)^ ■ l{\zn - c\ < 6)}. 
The corollary now follows from Theorem 12. 1[ □ 

Corollary 2.5. Under Assumptions AI-A4, for any < 6' < 6, 

ruA\zn-c\>6'] = O(n^i). 
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Proof. It follows from Chebyshev's inequality and Theorem 12.11 that 
FuA\^n-c\I[\zn-c\<6]>6'/2] < AEn^{\zn-c\^I[\zn-c\ <6]}/{6'f = 0{n-^), 
and that 

Fn„[\zn - c\ > 6] < EnA\zn-c\I[\zn-c\> 6]}/6 = 0{n-'), 
from which the corollary follows. □ 

3. The distance between n„ and its unit translation 

A key step in the argument leading to our approximation is to establish that 
the equilibrium distribution n„ of is sufficiently smooth. In order to do so, 
we first need to prove an auxiliary result, showing that, if the process Z„ starts 
near enough to nc, then it remains close to nc with high probability over any 
finite time interval. This is the substance of the following lemma. 

Lemma 3.1. Under Assumptions AI-A4, for any < t] < 6, there exists a 
constant Ku ^^ < 00 such that 

P[ sup \Znit) — nc\ > nrj \ Z„(0) = i] < n~^Ku^^, 
te[o,u] 

uniformly in \i — nc\ < nr]e~^^^ /2, where Ki := H-F'H^. 

Proof. It follows directly from Assumption A2 (a) that h defined by h{j) = j 
satisfies condition fll.l2p . Fix Z.„(0) = i, and define 

(3.1) := inf{t>0: \Zn(t) — nc\ > nrj}. 
Then it follows from Theorem 11.21 that 

Mnit) := ZnitAT^)-i- / nF{znis))ds 

Jo 

is a martingale with expectation 0, and with expected quadratic variation no 
larger than 

(3.2) nt J2 

jez\{o} 

at time t (see Hamza and Klebaner (1995, Corollary 3)); here, as earlier, Zn '■ = 
n~^Zn. Hence we have 

\zn{tAT^)-c\ < -\ sup \Mn{s)\ + \i - nc\\ + [ \F{zn{s))\ds, 
^ [selo,u] J Jo 

for any < t < U, and also, from Assumptions A1-A4, we have 

\F{z)\ = \F{z)-F{c)\ < sup \F'{y)\\z-c\. 

\y-c\<5 

Hence it follows that 

ptATri ptATri 

/ \F{zn{s))\ds < Ki \zn{s)-c\ds. 
Jo Jo 

Gronwall's inequality now implies that 

\zn{t A Tr,) - c\ < u'^l sup \Mn{s)\ + \i - nc\ > e^^\ 

s&[0,U] 
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for any < t < U, and so, for \i — nc\ < nrje^^'"^'^ /2, 

(3.3) sup \zn{tAT^)-c\ < ri/2 + n~^ sup \Mnis)\e^^^^ . 

t&[o,u] se[o,u] 

We have thus shown that 
(3.4) 

P[ sup \zn{t)-c\ > 7] I Z„(0) = z} < P[ sup \Mn{s) \ > ne-^^^7]/2 I Z„(0) = i]. 

te[o,c/] se[o,;7] 

But by Kolomogorov's inequahty, from (13. 2p . we have 
(3.5) 

P[ sup \Mn{s)\> ne-^^^'q/2\ Zn{Q) =i] < An-^'q-'^e^^'^U ^ fcj{l+ri), 

completing the proof. □ 
We can now prove the main theorem of this section. 

Theorem 3.2. Under Assumptions AI-A4, there exists a constant K > such 
that 

where n„ * 61 denotes the equilibrium distribution n„ of Zn, translated by 1. 

Proof. Because we have little a priori information about n„, we fix any [/ > 0, 
and use the stationarity of n„ to give the inequality 

ciTy{n„,n„ * 5i} 

(3.6) < ^n„(z)dTy{/:(^n(f/)|^n(o) = z),/:(z„(f/) + i|z„(o) = z)}, 

By Corollary 12.51 we thus have, for any 5' < 6, 

(3.7) dTv{Iln,Iln*Si} < Di„((50+O(n-i), 
where 

DmiS') := Yl n„(2) drvmZ^iU) \ Z„(0) = t), C{Z^iU)+l \ Z„(0) = t)}. 

i: \i—nc\<5' 

This alters our problem to one of finding a bound of similar form, but now 
involving the transition probabilities of the chain Zn over a finite time U, and 
started in a fixed state i which is relatively close to no. 

We now use the fact that the upward jumps of length 1 occur at least as 
fast as a Poisson process of rate A'', something that will be used to derive the 
smoothness that we require. We realize the chain Z^ with Z„(0) = i in the 
form Nn + Xn, for the bivariate chain {Nn,Xn) having transition rates 

(/, m) ^ (/ + 1, m) at rate nX^ 

(/,m) ^ (/,m + 1) at rate n Ai(^^) - A° 

(/, m) — > (Z, m + j) at rate nXj ^^-^ j , for any j G Z, j 7^ 0, 1, 
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and starting at (O,?). This allows us to deduce that 
dTv{C,{Z^{U) I Z„(0) = i),C{Z,,{U) + 1 I ZrM = i)} 

= \nZn{U) = k I Z„(0) =t)- F{Z^{U) = k-l\ Z„(0) 



IE 



2 



5^ P(iV„(f/) = mXniU) =k-l\ N^{U) = /, X„(0) = i) 

- 5^P(iV„([/) = / - l)P(X„(f/) = A; - Z I K,{U) = I - 1,X„(0) = 
i>i 

^ ^ E E = - nNniu) = I - i)\f[;Ak - o 
(3-8) + ^ E E = ^ - - - /.-m(^ - 01, 

k&Z 1>1 

where 

(3.9) /,;^(m) := P(X„([/) = m | N^{U) = /,X„(0) = z). 
Since, from Barbour, Hoist and Janson (1992, Theorem l.C), 

(3.10) |]|P(iV„(f/) = /)-P(iV„,(f/)=/-l)| < -^2= = 0(-^ 



the first term in (13.81) is bounded by 1 / {V nX^U} , yielding a contribution of 
the same size to Din{S') in (13. 7p . and it remains only to control the differences 
between the conditional probabilities fj^i{m) and fjLuijn). 

To make the comparison between /['ii'm) and flLn{m), we first condition 
on the whole Poisson paths of Nn leading to the events {Nn{U) = 1} and 
{Nn{U) = / — 1}, respectively, chosen to be suitably matched; we write 





I[o,uY 


. . dsi-i ds* 






P(X„([/) = 


--m\N^,[0,U] = u\- -si,.. 


,Sz_i,s*),X„(0) = z) 




I[o,uY 


.dsi-ids* 




(3.11) 


nxniu) = 


m 1 iV„[0,t/] =z/'-^(- 


. . ,sz_i),X„(0) = i). 


where 









V"" {U\tx, . . . ,tr) := ^ ]l[o,«](tj)) 

2=1 

and y [0, u\ is used to denote (l'(s), < s < u). Fixing si, S2, . . . , let Pj 
denote the distribution of conditional on A''„[0, \J\ = z/'(- ; Si, . . . , s*) 
and X„(0) = i, and let Pj denote that conditional on A^ri[0, U] = z/'"^(- ; si, . . . , Si_i) 
and X„(0) = z; let ps*{u,x) denote the Radon-Nikodym derivative ciPj,, 
evaulated at the path x[0,u]. Then 

P,,,4X„,(f/) = m] = / p,.(f/,x)rfP,(x[0,f/]), 

J{x[0,U]: x{U)=m} 
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and hence 
(3.12) 

¥,AXn{U)=m]-¥,[Xn{U)=m] = J l{^y{x{U)){ps4U,x)-l} dF,{x[0,U]). 
Thus 



1 / ds, ...dSi_,ds*y^E,{l{^y{Xr,{U))\ps4U,X,,)-l\} 

(3.13)< 4/ dsi ...dsi_ids*E.,{[l-ps4U,X^)] + }. 

J[o,uY 

To evaluate the expectation, note that ps*{u, Xn), m > 0, is a Pj-martingale 
with expectation 1. Now, if the path x[0, U] has r jumps at times ti < ■ ■ ■ < tr, 
writing 

y{v) := x{v) + u^-^{v ;si,...,si_i), yk := 2/(4), jk ■= Vk - Vk-u 
we have 

{1 if M < s*; 

exp (-n /; {X{y{v) + n"^) - X{y{v))} dv) 
U{k: s*<t,<u} {hkiVk-i + n~^)/Xj^{yk-i)^ if n > s*, 

where Aj(-) = Aj(-) if j 7^ 1 and Ai(-) = Ai(-) — A°, and where A(-) := 
J2j£Z\{o} ^j(')- Thus, in particular, ps*{u,x) is absolutely continuous except 
for jumps at the times tk- Then also, from Assumptions A3 (a) and A4, 



uniformly in |?/ — c| < 5, for each j G J. Hence it follows that, if we define the 
stopping times 

Ts := mi{u > : |X„(n) + i/~^{u ; Si, . . . , S;_i) — nc\ > n5}] 

(3.14) := mi{u>0: ps4u,Xn)>2}, 

then the expected quadratic variation of the martingale ps* {u, X„) up to the 
time mm{U, ts, 0} is at most 

(3.15) AU (i^) ncj{l + 6) =: n-^K{6,e)U, 

jez\{o} \ / 

where K{6,e) < 00 by Assumption A2 (a). 

Clearly, from fl3.15p and from Kolmogorov's inequality, 

Fi[(j)<mm{U,Ts}] < K{6,e)U/n. 

Hence, again from fl3.15p . 

Ea[l-P..(f/,X„)] + } < n-'/WK{6,e)U + n-'K{S,e)U + F,[ts < U]. 
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Substituting this into (13.131) . it follows that 



i>i 



< 2|n"i/ V^l*^. £)U + n-^K{6, e)U 

+ P[ sup \Zn{u) - nc\ > n6 \ Z„(0) = i]}. 

0<u<U J 

But now, for all i such that \i — nc\ < n6' = n(5e~^^^^/2, the latter probability 
is of order 0{n~^), by Lemma [3.11 and hence the final term in fl3.8p is also of 
order 0(n~^/^), as required. □ 

As a consequence of this theorem, we have the following corollary. 

Corollary 3.3. Under Assumptions AI-A4, for any bounded function f, 

EuAVfiZn)} = 0(^11/11). 

Proof. Immediate, because 

|En„{v/(^n)}| < 2||/||dry(n„,n„*5i). 

□ 

4. Translated Poisson approximation to the equilibrium 

distribution 

We are now able to prove our main theorem. The centred equilibrium distri- 
bution of Zn is n„ := n„ * 6^inc\, and we approximate it by a centred Poisson 
distribution with similar variance. 

Theorem 4.1. Under Assumptions A1-A5, 

dTv{^o{nv,),Iin) = 0(n-°/2), 

where := a'^{c) / {-2F\c)} . 

Proof. We follow the recipe outlined in Section [TTTl From (11. 9p . we principally 
need to show that 

sup|E{w \j g{W +l)-Wg{W) + {v)g{W)}\ = 0(n^°/2), 

for W := Zn — [nc\, v := nVc and E := En„. So, for any g G Gnv^, write 
g{i) := g{i — \nc\), and set 



0, if i < [ncj — [nvc\ ; 

Note that, for j > 1, by Assumption A2 (a), 

j 

n\j{i/n)\h{i + j) — h{i)\ < njcj||^|| + Cj\i — \nc\ \ \g{i + j — k — [nc\ 



k=l 

j 

< njcj\\g\\ +jcj sup \lg{l) \ + ^ |j - k\\\g\ 
' k=i 
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and that a similar bound, with \j\ replacing j, is valid for j < —1. From the 
definition of Qnvc fH-^p and fll.7p and from Assumption A2 (a), it thus follows 
that {\An\hn,g) is a bounded function, and hence that the function h„ g satisfies 
condition (11.121) : furthermore, since \hn,g{i)\ < |i— [ncj + [nwj I, in view of (11.71) . 
hn_g is integrable with respect to n„, because of Theorem 12.11 Hence it satisfies 
the conditions of Theorem 11.21 from which we deduce, as in (II. ip . that 

En„(A/in,3)(^n) = 0. 

Applying Lemma since hn^g has bounded differences in view of (II. 7p . it 
follows that 

= En„|^a2(^) v^(^n)+nF(^)^(ZO+K(^,^n)} 
= -F'(c)En„ {nvc V 9{Zn) - {Zn - [nc\)g{Zn) + {nVc)g{Zn)} 

(4.1) +En„{K(^,Z.)+^n(^,^n)}, 

where En is as defined in (II. 3p . and 

KM ■■= -{a\^/n)-a\c))^g{^) 

+ {n{F{i/n) - F{c)) - F'{c){z - [nc\)}g{z) + F'{c){nv,)g{i). 

The terms involving E'^{g,i) can be bounded, using (ll.7p . as follows. First, 
using Assumptions A2 (a) and A4, 

l\a'{z/n)-a\c)\\vm\ 

< II (cr^)'||5|z — nc|/[|z — nc| < n5] 

(4.2) +^( E fcj{l + Wn-c\) + a\c)y[\i-nc\>n6]- 

j6Z\{0} 

and then, under Assumptions A2 (a) and A5, 

\n{F{t/n) - F{c)) - F'{c){t - [nc\) + F'{c){nvc)\ \g{i)\ 
= n\F{t/n) - F(c) - {i/n - c)F'{c)\ \g{i)\ 

< ('^{i/n-cyi[\i/n-c\<6] sup \F"{z)\ 
(4.3) +nUl + \i/n-c\) ^ \j\cj + F' {c)\i/n - c\\l[\i - nc\ > 6]) 

iGZ\{0} ^ 

The contribution to (14. ip from En„{-E'^(^, Zn)} is thus of order 

En„{|z„ - C| + (1 + \Zn - c\)I[\Zn - c\ > 5] + \Zn - c\'^ I[\Zn " c| < 5]} 

(4.4) = 0{n-^'^), 

by Theorem 12.11 and Corollaries 12.41 and 12.51 The first term in En{g,i) is also 
bounded in similar fashion: from Assumptions Al, A2 (a) and A4, 

^\F{z/n)\\vm\ 

(4.5) < ^{||F'||5|i-nc| + Cj\j\{l + \i-nc\)I[\i-nc\> 6]}. 



i6Z\{0} 
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giving a contribution to E,Yi^{En{g, Zn)} of the same order. The remaining 
terms, involving \j'^g, need to be treated more carefully. 

We examine the first of them in detail, with the treatment of the second 
being entirely similar. First, if either \i/n — c\ > 6 or j > ^/n, it is enough to 
use the expression in f ll.4p to give 

(4.6) \a,Cg,t)\ < j(j-l)IIV^II < j(j - 1)/Hc). 

For \i/n — c\ > 6, by Assumption A2 (a), this yields the estimate 



I[\i - nc\ > 5] 



(4.7) 



^aj{g,i)nXj{i/n) 
< J2^^^^^^^{l + \t/n-c\)I[\t-nc\>6], 



with corresponding contribution to En„{-E'n(fi', ^n)} being of order 0{n ^), by 
Theorem 12.11 and Corollary 12.51 Then, for j > ^/n and \i/n — c\ < 6, (14. 6 p yields 



^ aj{g,i)nXj{i/n) 



(4.8) 



< $:^^^^(l + 5) < 5^/+%n-/^(l + 5)/t;. 



making a contribution of order 0{n~°'/'^) to En„{-E'n(^, -^n)}, again using As- 
sumption A2 (a). In the remaining case, in which j < ^/n and |i/n — c| < 5, we 
use ( II. 5p . observing first that 

nsy'^ g{i + J - k + l)\j{i/n) 

(4.9) = ns^^g{i + j-k + l)\j{c) + n\7^g{i + j-k + l){Xj{i/n)-\,{c)), 

the latter expression being bounded by 



(4.10) \ns^^~g{z+j-k + l){\,{z/n)-\,{c))\ < -\\\'AU\^/ 



n — c\. 

The corresponding contribution to Ku„{En{g, Zn)} is thus at most 

J2 (jV6){A,(c)nsup |En„ v' K^n + l)\ + 2v;'\\X'^\\sEuJz^ - c\} 

3=2 ' 

< ^ j2+"c,{nsup |En„ v' 9{Zn + 01 + Li2v^'EuJzn - c|} 

i>2 ' 
(4.11) = n^'-''^/^0{n-n-'/^ + n-'/^) = ©(n^/^), 

where we have used Assumptions A2 (a) and A4, and then Corollaries 12.41 
and 13.31 and finally (11. 7p . 

Combining the bounds, and substituting them into (14. ID . it follows that 

|lEn„ {nvc V 9{Zn - [nc\) - {Zn - \nc\)g{Zn - \nc\) + {nVc)g{Zn - \nc\)] \ = 

uniformly in G Qnvc ■ Again from Corollary 13. 3^ we also have 

\nVcEn„ {V9{Zn - [nc\) - \jg{Zn - \nc\ + 1)} | = 0{n-^'^), 
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for any g G Qnvc- ^ thus follows from (11.91) that 

rfry(PSHe), fin) = O (n-"/2 ^ p^j^^ - nc < - [nv,\]) , 

and the latter probability is of order 0{n~^) by Corollary 12.51 This completes 
the proof. □ 

Example. Consider an immigration birth and death process with births 
occurring in groups of more than one individual at a time. The process has 
transition rates as in Section 11.11 with 

A_i(z) := dz, Ai(2:) := a + hq^z and \j{z) := hqjZ, j > 2, 

while Xj{z) := 0, j < —1. Here, b denotes the rate at which birth events 
occur, and a > represents the immigration rate. The quantity qj denotes the 
probability that j offspring are born at a birth event, so that J2j>i Qj = Ij 
write rrir := Ylj>ij^Qj '"'^'^ moment of this distribution. Then 

F{z) = a + z{hmi — d), and (j'^{z) = a + z{bm2 + d). 

Assumption Al is satisfied if c? > bnii, with c = a/{d — brrii) and F'{c) = 
— {d — hrrii). Assumption A2 (a) is satisfied with Cj = feg^ max{l, c}, j > 2, 
Ci = max{6gi,a + bqic}, and c_i = dmax{l,c}, provided that m2+a < oo for 
some < a < 1; for Assumption A2(b), simply take A" = a/2. The other 
assumptions are immediate. 

The quantity Vc appearing in Theorem 14.11 then comes out to be 

a(2d + 6(m2 — mi)) 
2{d-hmiY ' 

and the approximation to the equilibrium distribution of Zn — Ync\ is the centred 
Poisson distribution Po(nt;c), accurate in total variation to order 0(n""/^). 
Note that, if 6 = 0, then the process becomes a simple immigration death 
process, whose equilibrium distribution is precisely the Poisson distribution 
Po (na/(i) = Po (nc). In this special case, the approximation is in fact exact. 
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